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Abstract 



We analyze the scalar field sector of the Kazakov-Migdal model of induced 
QCD. We present a detailed description of the simplest one dimensional 
(d—1) model which supports the hypothesis of wide applicability of the mean- 
field approximation for the scalar fields and the existence of critical behaviour 
in the model when the scalar action is Gaussian. Despite the ocurrence of 
various non-trivial types of critical behaviour in the d—1 model as N — > oo, 
only the conventional large- N limit is relevant for its continuum limit. We 
also give a mean-field analysis of the N = 2 model in any d and show that a 
saddle point always exists in the region m 2 > m 2 rit {= d). In d — 1 it exhibits 
critical behaviour as m 2 — > m 2 rit . However when d>l there is no critical 
behaviour unless non-Gaussian terms are added to the scalar field action. 
We argue that similar behaviour should occur for any finite N thus providing 
a simple explanation of a recent result of D. Gross. We show that critical 
behaviour at d>l and m 2 > m 2 rit can be obtained by adding a logarithmic 
term to the scalar potential. This is equivalent to a local modification of the 
integration measure in the original Kazakov-Migdal model. Experience from 
previous studies of the Generalized Kontsevich Model implies that, unlike the 
inclusion of higher powers in the potential, this minor modification should 
not substantially alter the behaviour of the Gaussian model. 



1 Introduction 



The Kazakov-Migdal model QXJ| for induced QCD (to be further referred to as 
the KMM) has recently attracted much attention ||-||. It is interesting as 
a tractable example of a matrix model where the the integration over angular 
variables plays a non-trivial role and as such, it is the natural next step in the 
theoretical investigation of matrix models. It is also a lattice gauge theory 
and there is the possibility that, provided the appropriate scaling limit exists, 
its large wave-length behaviour is described by Yang-Mills theory. The KMM 
contains two kinds of fields, NxN Hermitean matrices $ which live on the 
sites of a d- dimensional lattice and NxN special unitary matrices U which 
live on links. At large distances, the latter are associated with the gauge 
fields of Yang-Mills theory, while the former play an essential role at short 
distances and presumably disappear from the spectrum in the continuum 
limit. 

The most direct approach to finding a continuum limit of the KMM is 
to assume that N is large and to use the mean-field approximation for the 
$-fields. It is natural to assume that the mean field is constant in space- 
time. The first results in the framework of this program were obtained in 
[0]. There an equation for the density of eigenvalues, p(<f>), of the master- 
field $ was deduced from a highly non-trivial saddle point equation with 
the help of an especially simple subset of the Schwinger-Dyson equations. 
Using a power-like ansatz, p(<f>) ~ 0°, a solution of this equation was found. 
However, in |7j an exact solution of the saddle point equations was found 
for the case when the action for the scalar fields is quadratic. In that case 
the eigenvalues had a semi-circle distribution vrp(0) ~ \J1p — p 2 <p 2 . It was 
further argued that for d > 1 there is no critical behaviour (p — > oo) if the 
critical point is approached from the strong coupling phase (m 2 > m 2 rit ). 
One must rather approach the critical point from the weak coupling phase 
(m 2 < m 2 rit ). In order to make the effective potential for the scalar field 
stable in this situation one needs to introduce higher order (at least quartic) 
terms in the bare action. The resulting critical behaviour will generically 
be first order though second order critical points may be present for some 
special values of the parameters. The N = 2 version of the KMM with 
quartic terms in the potential has been examined in the region m 2 <m 2 rit 
in [[J where such critical behaviour was indeed found. Moreover the results 
of computer simulations demonstrated a nice agreement between the N=2 
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KMM and the predictions of the mean field analysis which has no a priori 
reason to be reliable for N = 2. 

In this paper we shall present a more extensive analysis of the easily 
solvable versions of the KMM. We believe that this is necessary before aban- 
doning the simplest version of the model with a quadratic potential for the 
scalar field or accepting the suggestion (implied in || and 0) that there 



may be second order critical behaviour in the instability region m 2 <m 2 
(We shall call the model with such a quadratic potential a "Gaussian" model 
despite the fact that the potential is only quadratic in the scalar ($) fields. 
Its interaction with the gauge fields ([/-variables) is highly non-linear.) 

We begin by analyzing the exactly solvable 1 -dimensional KMM in detail. 
There we observe some amusing types of critical behaviour which occur when 
iV — > oo but which do not survive in the thermodynamic limit (when the size 
of the system increases to infinity). We shall also use the explicit results in 
d=l to demonstrate the applicability of the mean field approximation both 
for these different kinds of critical behaviour as well as for the simplified 
N = 2 version of the KMM. The mean-field results for the continuum limit 
of the N = 2 model are then generalized to any d > 1. In this case the 
critical behaviour disappears in the region m 2 > m 2 vit = d (and survives 
only for m 2 < m 2 rit ). Thus it seems necessary to deal with the upside-down 
harmonic oscillator potential and possibly to introduce higher order terms 
in the scalar potential. We then argue that a similar conclusion applies to 
all N>2. This provides us with a more transparent explanation of the result 
of ref. |7| concerning the absence of criticality in the Gaussian model when 
d>l. The origin of this phenomenon is just a logarithmic increase of the 
effective potential at infinity due to the Van-der-Monde determinants. It is 
possible to compensate for this growth by introducing logarithmic terms into 
the bare potential. These can be interpreted as a change of the measure of 
the integration over matrices. Changes of the measure of this kind have been 
investigated in the somewhat simpler context of the generalized Kontsevich 
Model |J in which case it is known to preserve all of the nice features of 
the model with a quadratic potential []T0| | . This may also be the case for the 
KMM though further investigation is required. 

Our analysis of the various types of critical behaviour (even those irrele- 
vant for the continuum limit but, rather, associated with the large N limit 
in finite volume) indicates a wide applicability of the mean-field approxi- 
mation for the ^-fields. The result for the N = 2 model also supports the 
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suggestion that the critical behaviour for m 2 — > m 2 rit + (whenever it occurs) 
is associated with a mean field $ which is large, i.e.$ ^> 1. (This is not nec- 
essarily true if the critical behaviour is associated with the tuning of higher 
order terms in the bare potential for $.) Both of these features provide some 
justification for the assumptions which were made in |J in our analysis of 
the correlators of physical observables in the KMM. 

We begin in Section 2 with a review of the KMM. Sections 3 and 4 are 
devoted to a detailed discussion of the 1-dimensional KMM. In Section 5 
we consider the N = 2 model in any dimension and investigate its critical 
behaviour in the mean field approximation. 

2 KMM: Some Generalities and the Mean- 
field Approximation 

The KMM contains Hermitean matrices $(x) and special unitary matrices 
U(x, y) which are defined on the sites (denoted by x) and on the links (de- 
noted by < x, y >) of a d-dimensional lattice, respectively. The partition 
function is given by 

Z D = [d$[dU]exp I -^tr1/($(s)) + ]T tr$(x)U{x, y)${y)tf{x, y) ) (1) 

where [dU] is the invariant Haar measure for unitary matrices. For most of 
this Paper we shall consider the case where the action for the scalar fields 
is quadratic i.e. V A ($) = m 2 $ 2 . Higher order terms could be used to cut 
off the fluctuations of large $-fields near the critical point and would affect 
the detailed properties of the continuum limit. We shall argue that in d=l 
these additional terms are unnecessary. We shall also argue that in d>l 
at least some additional logarithmic terms (which could also be viewed as 
modifying the measure rather than the action in ([I])) are needed to produce 
the appropriate critical behaviour. 

If the integral over $ is taken first in eq.(^) it gives rise to an effective 
action for the gauge fields 

9 (U)- VMt (2) 
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where the sum is over all loops T and U[T] is the product of [/-matrices 
associated with links in T. Besides the conventional symmetry under gauge 
t r ansf or mat ions , 

$(x)^>V(x)$(x)V*(x) , U(x,y) -»• V(x)U(x,y)V^(y) (3) 

with V(x) G SU(N) there is also an invariance under multiplication of link 
operators by an element of the center of the gauge group || , 

U (x, y) -> uj(x, y)U(x, y) , (4) 

where, for the U(N) group u(x,y) is a unimodular complex number and 
for the SU(N) group it is an element of Z N , i.e. (u(x,y)) N = 1. This 
symmetry has important implications for the properties of observables and 
their correlators |J. 

In the "naive continuum limit" this effective theory resembles Yang-Mills 
theory provided £) = 4 0. However, this limit is too naive in the sense 
that, unlike the conventional Wilson formulation of lattice QCD which has 
a single, or finite number of terms like 

^tvU(V) (5) 

in the action, the effective action (0) has the special property that it is 
not sharply peaked in the vicinity of the trivial configuration U = 
This leads to significant differences between the procedure for taking the 
continuum limit in this gauge theory and in conventional lattice QCD. 

Remarkably it is possible to do the gauge field integral first in (p. The 
result is an effective action for the $-fields, 

S cfr ($)=m 2 5>$ 2 (x)- J2 log/($(x),$(sO) (6) 

x <x,y> 

where 

denotes the integral over the [/-matrix on a given link and we denote by 
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Vn = Vol(Z7(iV)) the volume^ of the unitary group (in the Haar measure 
[dU\): 

N 

V^=(27r)^/n*! (8) 

k= 

The integral (0) depends only on the eigenvalues of the ^-fields, denoted here 
by 4>i and 

a(0) = n(& - (9) 

i<j 

is the Van-der-Monde determinant. This explicit formula for \E r ) has 
been known for a long time [11]. In the large N limit the resulting effective 
theory for the $-field can be analyzed in the mean field approximation. 



4 A simple way to derive this formula is provided by the theory of random matrices. 
The Gaussian integral over Hermitean matrices H , 

JV 



Jn = JdHe~ UH2 =f[ J dH u e^l[J dH l3 dH l]e - 2 \ H 



A 2 = n N ft 



i— 1 i>j ' 



can alternatively be expressed by diagonalizing H, integrating over orbits of diagonal H 
and the orthogonal polynomial method to compute the remaining integral as 

„ N 

{ \ — T JV 9 

J N = Vol[U(N)/U(l) N ] / Y[dhi& 2 (h)e~^*=i h * = 

i= 1 

„ JV N-l 

Vol[U(N)/U(l) N ]N\ / Y[dhie-^ h i J] **i fa) 

J i=l j=0 

The last product on the right-hand-side contains the norms of orthogonal Hermite poly- 
nomials normalized so that 

H k (h) == ±e h ''*(h ^) k e- h2 / 2 = h k + ... 

Then ||-H/c|| 2 = \pHk\j2 k . The JV! factor in front of the integral comes from the 
permutations of the eigenvalues hi while the volume factor arises from integration over 
the flag manifold U(N)/U(l) Ar of unitary matrices U(N) modulo the Cartan (diagonal) 
elements of U(N) which commute with the matrix H once it is diagonalized. Comparison 
of these two expressions for gives (||). 
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The main idea of ref . was to analyze the classical equations of motion 
for the effective action 5 , e if($), 

m®{x) = \ £ f «^"t) log (10) 
together with some simple Ward identities for /($, 

tr (^) J ( $ '^) = tr ^ fe/ ( $ '^)- ( U ) 

It is important to note that this is not a complete set of Ward identities. 
They do not define /(<!>, \l/) unambiguously. For example 

J ($,*) = e tr ** (12) 

is also a solution of (THD) . A complete set of identities would contain analogues 

of these equations with operators like tr (^$) \ (where Z^O) on the left 

hand side). There may also be independent identities with mixed $ and \I> 
derivatives. Unfortunately most of them (with the exception of ([11])) have a 
complicated form. 

Since the reduced set of identities (|TT|) is incomplete, it is not clear 
whether arbitrary solutions to the equations [§] 

K(\) = P f°*l w a -^)-^^)+^) = 

1 ; J-oo 2m A - ^V'{v) - ^j-H{y) - iirp(u) 

+°° dv ^piy) 

— arctan = — -3— ; — — — (13) 

with 

W> s pfH®±, (14) 

derived from them will have anything to do with the actual KMM.0 While the 
results of Migdal and Gross might still be relevant to the KMM for some less 

5 In any case it is unclear how to apply the same trick to estimate quantities other than 
the eigenvalue distribution. To study correlation functions in the KMM one must use some 
more substantial information about The feature ( pT| ) is insufficient (see Q). 
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obvious reason this reason needs to be clearly identified. It was shown in M 
that ( |I3| ) posses a solution for p(/i, d) in the form of a semi-circle distribution 



np(ji, d) = yj2ti(d) - ^(d)^ (15) 



which is usually characteristic of matrix models with purely quadratic po- 
tentials. This is not surprising since the information input in the derivation 
of (p~3|) does not distinguish between the integral J($, \P) in (|7|) and the triv- 
ial quadratic potential 7 (^ ) ,^) in (0)- One reason why ( jl3|) may still be 
relevant is that for large 0, (see (|59[)), /($, \l/) is not very different from 



was derived in by substituting the semi-circle ansatz (|T5| ) into ([13]) .Q This 



expression shows that while for d = 1 // goes like ym 2 — m 2 T - lt in the vicinity 
of the critical point this is no longer true for g? > 1 in which case the branch of 
fi which is positive in the region m 2 >m 2 rit = d does not vanish when m 2 =d. 
We shall see in Section 5 below that this result has a very simple explanation 
and we suggest a way to restore the critical behaviour for d>l. 

3 Partition Functions for Various Models in 
One Dimension 

In one dimension the KMM can be solved using elementary techniques. It 
essentially reduces to the well known problem of a matrix-valued harmonic 
oscillator. This problem has recently been reviewed in the context of matrix 
models in the recent papers [0, j| as well as in || and 0. There are 
as many as four slightly different Gaussian matrix models in one dimension. 
They can be defined either on an open or on a closed lattice (chain), and they 
can either contain only hermitean matrices as in conventional matrix models 
or they can couple to unitary matrices to form a gauge invariant model. 

We shall now list the four models and give the expression for the partition 
function for each. The relevant parameter in this Gaussian problem is q± = 

6 Here we shall use our notation, for example our m 2 is two times smaller then in M 



I ($,¥). 

The explicit expression 




(16) 
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m 2 ± -y/m 4 — 1 (rather than the "bare mass" m). In the closed chain models 
logg is additive in the length L of the chain. We thus define qi = q L ■ In the 
discussion below we use the branch q_. 



3.1 Gaussian Multimatrix Model on an Open Chain 

This model contains a chain of coupled matrices with free boundary condi- 
tions. The partition function was computed in |5|: 

r L L-l //-i _ rl 2\ n L\ N2 / 2 

Z A ^jU d^e- 2 ^ n = (17) 

^ x=i s=i \ 1 y / 

3.2 Gaussian KMM on an Open Chain 

The partition function for KMM on an open chain is essentially the same as 
the partition function for the previous case. 

Z B = J f[d® x e~ mW *]l[dU x>x+1 ^ 

x=l x=l 



//i 2\ L\ N I 2 

= [Vol(U(N))] L Z A = V£ [ \: q q2 [ q +2 ) (18) 

The reason for this is that the gauge transformations (§) can be used to 
substitute unit matrices / for all of the U XtX+ i variables on the chain. The 
resulting integral is identical to that of the Gaussian multimatrix model with 
the exception of the factor Vn which is the volume of the group U(N) with 
the Haar measure [dU]. 



3.3 Gaussian Multimatrix Model on a Closed Chain 

We now consider the multimatrix model on a closed chain with periodic 
boundary conditions = $i- The partition function for this model can 

be easily evaluated as 

^^.-s^^-fe (19) 
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3.4 Gaussian KMM on a Closed Chain 



For the Gaussian KMM on a closed chain the partition function is given by 



r L 

J$r_i_i=$l _ i 



X=l 



7V 2 /2 

= v *mlrw) (20) 

In this case the gauge transformations (||) can be used to eliminate (i.e. 
put equal to I) all the [/-matrices except for one. This last U can only be 
diagonalized so that U a b = e l9a 5 a b. Using this idea one can derive another 
expression for Z D M,M' 



N2 / 2 , n n ^ _ p i(e a -e b ) 2 



z D = {2,rv^ J ^— w J\{de a \[ 

\ L HL) •> a= i i<j 



(21) 



The right hand side of (|20D can be considered as the result of performing the 
integration in (Elf). 

3.5 Comments and Applications 

We now make several observations which are important for our further anal- 
ysis of the d=l model as well as for its generalization to higher dimensions. 
First of all the factor 

?r 12 = <i ln2 > 2 (22) 

is common to all four partition functions. The main difference between them 
is their behaviour at the critical point g crit = 1 (i.e. m^ it = d = 1). Note that 
both the functions Za and Zb for the open chain models are not singular at 
this point. On the other hand the partition functions for the closed chains 
do possess a singularity as q — > 1. The nature of this singularity is, however, 
quite different in the two models. In the conventional multimatrix model 

Z c ~ (1 - qY N ' 2 as q -> 1 (23) 
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while in the KMM 



Z D ~(l-q)- N asg^l (24) 

The presence of this singularity reflects the occurence of zero modes in the 
integral as m 2 — > d—1. In the model C, any constant (x-independent) matrix 
= $ is a zero mode in this limit. This results in the power N 2 in 
(|23|). Things are however different for the KMM D. Integration over the 
[/-matrices can, and indeed does, eliminate some of the zero modes which 
occur at U — I. (This was the essence of our argument in |J against the 
"naive continuum limit" for the [/-variables which would imply that the 
vicinity of U = I gives the main contribution to the integrals over U.) The 
key observation is that in fact not all of the zero modes disappear. The zero 
modes in Cartan subalgebra survive. Indeed if the [/'s can be diagonalized by 
a gauge transformation (^) (as we have seen is the case in the 1-dimensional 
model) then the diagonal elements of the matrices $ do not feel f/'s. In fact 

tr$[/W f -> J2 ^ab^bae i{da ' eb) (25) 

a.b 

so that all coupling between the diagonal elements of $ and \l/ and the #'s 
disappears. This fact, combined with the absence of coupling between the 
diagonal and the off-diagonal elements of $ and \1/ in ( |23| ) ensures that the 
zero modes $ aa (x)=$ aa = const cannot be eliminated by the [/-integration. 

We can understand how the off-diagonal zero modes are eliminated by 
looking at eq. fl21~|) . Note that the integrand is not singular at q^ = 1. The 
singularity coming from the action is completely canceled by the zeros of the 
measure. (This is a very unusual situation which is specific to the KMM. It 
can be understood as being due to the very slow i.e. logarithmic increase of 
the action near its minimum. This can in turn be traced back to the fact 
that contours of arbitrarily large length contribute to the effective action 
for [/-variables. See || for details.) The net result is that the order of the 
singularity in (|24"D is just iV = rank U(N). The generalization of this result 
to d>l is discussed in Section 4. 

We now briefly discuss the implications of this singularity. According to 
M the fact that the singularity is not of the power iV 2 in the KMM (as it 
was in the model C) implies that there is no "naive continuum limit" for the 
[/-variables (It is still possible that a "less naive" continuum limit may exist.) 
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Note however that the singularity still does survive. In fact it is of order N 
which, although much less than N 2 , is still much greater than 1. In the next 
section we shall argue that this fact implies that non-naive large N limits may 
exist in the $-sector even for finite volume. The resulting models are of no 
interest when taking the continuum limit since they do not survive in the limit 
of large volume. They are however interesting examples of critical behaviour 
and they should thus be kept in mind in future investigations. We shall also 
see that these limits are consistent with the mean-field approximation for $ 
and they thus provide us with additional material to study its features. 

4 Different Types of Critical Behaviour for 
the Gaussian KMM in One Dimension 

Our goal in this section is to discuss the various kinds of critical behaviour 
which are possible in d= 1 and to study the applicability of mean-field theory 
in each case. Since we have exact expressions for the various partition func- 
tions we can easily test the validity of the mean-field approximation. The 
main effect which results when mean-field theory is exact is the factorization 
property for correlators such as 

< (tr$ 2 (si) ■ • ■ tr$ 2 (x fc )) > =(<tr$ 2 >)\ (26) 
These particular correlators are especially suited for our purposes since 

«(^X>*W)\> =l(-^)"z. (27) 

Thus if the factorization property (|26| ) is exact we should find, in the large 
N limit, that 




Note that we have introduced a normalization factor M in the definition 
of these operators. The conventional choice (corresponding to the "naive 
continuum (or large N) limit" which was examined in f7|) is M nc i = 1/N 2 . 
However our experience in |J indicates that one should be very cautious 
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when fixing the value of Af. The correct choice of M can (and does) depend 
on precisely how the limit is taken. 

To illustrate the various possibilities for critical behaviour let us proceed 
directly to the most interesting case - the KMM on a closed chain (model 
D from the previous section). The partition function Z is given in eq. (]20|) . 



We begin by evaluating the correlators (|28|). Since Z depends on qL=q L and 
we wish to differentiate it with respect to m 2 we must use the fact that 

d\ogqL dlogq 1 2q 



L dm 2 dm 2 \Jm 4 — 1 1 — q 2 



(29) 



Then 
and 




N 



<<mi e» = M T ^ T2 {^ + Et^J (») 



n=l 




71=1 




We are now ready to discuss the different limits. We are concerned with 
the limit N — > oo, q — > 1 and possibly L — >• oo. Let us define 6^=1 — ^=1 — 
(1 — e) L and e = ei = 1 — q. We shall always be concerned with the limit 
e — » so that e is assumed to be small. We shall also assume that ei<^l so 
that eL<l. 

4.1 Limit (a): iVe L > 1 

Let us first consider the case when iVe^ ^> 1. This is the naive large iV limit 
(considered in 0). In this limit we pick up only the terms with the highest 
powers of N. If we choose the normalization M = 1/N 2 these are the terms 
which are finite in the limit iV — > oo. From the above formulae we find in 
this limit 

(M \ k ( AfN 2 \ k 
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(<A/tr$ 2 >) fe (l + 0{l/N)), (32) 

We thus see that the factorization property is valid in this limit. Further- 
more the particular value of <C 7Vtr$ 2 ^>= \/2\Jm 4: — 1 as well as the other 
averages (see below) are consistent with the prediction of the semicircular 
distribution 7rp(0) = \J1e — e 2 </> 2 0. 

4.2 The limit (b): N -> oo, iVe L <C 1 

The limit A" — > oo with A^e^ <C 1 can be quite different than the previous 
one since terms in ([H]) which have less powers of iV can be more singular as 
€l — > 0. In fact since 



N nn n N n N 

Er^ = E-(l + 0{ne L )) = - (1 + 0(Ne L )) , (33) 



and 



we have: 



N n 2 n n N 

£ra^ (1 + 0(JV£j) (34) 



Af (N 2 N N 
< A/tr$ 2 > = — — + - (1 + 0{Ne L )) 
e L \ 2 e L 



^ (1 + 0(Ne L )) (35) 



and 



^ (1 + 0(AT ei )) + — ^— + | (1 + 0(iVe L ))J | = (36) 
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(f)>+^™>) 

We see that in this limit the results are somewhat different. Firstly in or- 
der for a reasonable limit to exist M should be equal to 1/N (rather than 
1/iV 2 ). Secondly the averages behave as e L 2k (rather than e L k ). But despite 
these differences the factorization property still remains valid for fc C iV. 
Note however that the relevant master field is clearly different from that of 
ref . @ and there is no reason to believe that it is described by a semicircular 
distribution. Our analysis so far which is based on the calculation of the 
specific correlators <C (X)tr$ 2 ) ^> is not sufficient to completely rule out 
this possibility. 

To demonstrate the failure of the semicircular distribution, one should 
evaluate the more complicated averages <C A4tr($ 2fe ) ^> which cannot be 
directly extracted from the above partition functions. If the distribution of 
eigenvalues obeys the semicircle law these averages should be equal to 

< A4tr$ 2fe > semi circular= T^U ( 37 ) 

Any deviation of these averages from this formula indicates that the distri- 
bution is not semicircular. We compute these averages below. 

Before we proceed to the calculation note that the normalization factor 
Afc is equal to 

N k = ^(NAf) k . (38) 

This will be true for all the types of limits which we consider. The reason 
for this is that every trace should be accompanied by a factor of 1/N if it 
is to be substituted by an integral with an eigenvalue-distribution function 
jjtr... — > J d(f>p((f))... normalized so that / d<fip(<f)) = 1. The factor (NAf) k can 
be absorbed into the normalization of the $-field in which case the role of the 
effective Plank constant in the integral would be played by (NJ\f)/N = M 
(the factor of N in denominator is associated with traces in the action). 
Applicability of the saddle point approximation requires M to be small. 

To simplify our evaluation of <C A/j-tr$ ^> we shall begin with the case 
L=l. It is straightforward to generalize to the case for general L. Up to terms 
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which are negligible since they are of higher orders in 1/L and 1/N, it can 
be obtained by substituting qi for q. The integral of interest || 

J d9 a A 2 (e i6 °) J d^expi-YsMim 2 ~ cos9 ab )) tr$ 2fc (39) 

a, b 

can be greatly simplified in the vicinity of the critical point (This integral 
should, of course, be divided by a similar integral with k=0 in order to give 
<C A/fctr$ 2fc 3>). Indeed, for k=0 this integral is just equal to (0), and the 
integrand is non-singular at q—1. At the limiting point q—1 the integral 
is simply equal to unity It follows that the measure of integration over d9 
in (|39"D becomes trivial at the critical point. Thus in order to compute our 
correlators we need only evaluate the contributions of propagators 

< <Pab<Pcd > = (2(m 2 - cos e ab )Y l 5 ad 5 bc (40) 

Since we are interested in the most singular terms at small e = 1 — q ~ 
^2{m 2 — 1), we can take all the differences 9 ab = 9 a — 9 b in the propagators 
to be small, so that 

<0a60 C£i > ^(9 2 ab + 2(m 2 -l))- 1 5 ad 5 bc (41) 

The propagator for the diagonal components of $ is independent of 9 and 
equals 1/e 2 . The contribution of any off-diagonal component 

~ / d9 ab (9 2 ab + e 2 )" 1 / J d9 ab = (7r/2e)(7r)- 1 = l/2e (42) 

It remains to compute the number of diagonal and off-diagonal components 
arising from tr$ 2fe . Their contributions are easily distinguished since those 
of the diagonal components give rise to higher powers of 1/e but to lower 
powers of N. There are thus iV diagonal components and iV 2 (l + 0(1/N)) 
off-diagonal ones. Let us first consider the limit (a) where we retain only the 
terms of order iV 2 . Now, since k iV we can assume that integrations over 
different off-diagonal elements are independent. Using a combinatorial argu- 
ment (see below, this gives the contribution of the off-diagonal components 
to the correlator as: 

Af k s° k N k+1 /(2e) k (43) 
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with s° = (2k)\/k\(k +1)!. With the proper choice of normalization factor, 

TV = 1/N 2 (44) 

thus 

A4 = 1/N k+1 (45) 

this coincides with ( j3^) , and confirms the result of |7[] for the limit (a) in one 
dimension. 

In the limit (b) we instead take into account the contribution of only the 
diagonal elements of $. In contrast to the off-diagonal case, the contribution 
to tr$ 2fe is given by the N times the average of one component to the power 
2k. This gives a factor of N (in accordance with the proper normalization 
M = 1/N in this case) times 

d0(0) 2fe e- £ ^ 2/2 ~ (2k - l)!!e- 2fc = s k k e- 2k (46) 

In this limit we clearly get a distribution function of the form 

M = 4re-^ (47) 



rather than the semicircular distribution. Note that this p(4>) is non- vanishing 
on the entire axis. 



4.3 Limit (c): Ne L ~ 1 

We shall now consider the limit where N is large, e is small but Ne is of 
order 1. This is a kind of a "double scaling" limit somewhat different from 
the type usually encountered in matrix models of gravity where N and L are 
correlated and iV can be interpreted as a new continuum variable giving rise 
to a new space-time dimension d — > d + 1. The formulae for the correlators 
in this limit are 

<A4tr$ 2fc > = {NH) k j2 s K ^ N/ k 2 S \ l + 0(1/N)). (48) 

3=0 6 3 

In the limits (a) and (b) the first (j = 0) and the last (j = k) terms in 
this sum are dominating respectfully. The coefficients s% can be defined 
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from the following combinatorial problem. Take tr0 2fc = J2{a t } 4>a 1 a 2 --- ( t>a 2k a 1 
and consider all the pair contractions of $'s (Wick rule), each contraction 
being < (f) ab (f) cd >= 6 ad 5 bc (a + (35 ab ), where a = l/2e and (3 = 1/e 2 . The 
is by definition the leading term (with the highest possible power of N) in 
the coefficient in front of a k ~^ ft . Let us draw a picture (Fig. la) where dots 
correspond to the ^-matrices. Every solid-line contraction (Fig. lb) of dots, 
gives rise to a factor of a and makes an identification a = d,b = c; while 
every wavy-line contraction (Fig.lc) contributes a factor of /3 and identifies 
a = b = c = d. One can easily check that any intersection of two solid 
lines "eats up" a factor of iV 2 , while any intersection of a solid line and 
a wavy line eliminates one N. Thus the leading-iV contribution arises in 
the "planar" limit, when the only lines of contraction which are allowed to 
intersect are wavy lines. Thus, estimation of is a clear combinatorial 
problem of enumeration of pair contractions of the 2k points on a circle by 
k—j solid and j wavy lines, such that only wavy lines are allowed to intersect. 
One can easily derive a recurrent relation: 



s&i^TT^-E E'4 j) stP; *<k. (49) 
« + 1 - % 3=0 l=j 

It is enough to note, that the first solid line separates the circle into two 
new circles of the lengths / and k — I. The factor arises because there 

are 2{k + 1) points to which the first solid line can be attached, while at 
the end there are as many as 2{k + 1 — i) points at the ends of all solid 
lines. The relations fl4"9]) should be supplemented by "initial conditions": 
= 1; = (2k — 1)!!, the first one is obvious, the second is just the 

number of all possible (intersections allowed) contractions of 2k dots by k 
wavy lines. This information is enough to find any s^. In particular, 

(0) _ k (2k-m _ m 

Sk ~ Z (Jfc+l)! "*!(* + !)! [M) 



K-i+j 



and 

-A[k] (51) 



,(i) = p i\k] _ (2fc)! 

' fc Sk (k-i)\i\ (k + 2)...(k + i) {k - i)\{k + i)HY 



where Pi[x] are certain polynomials of degree i-1: 



Po = ^, Pi = l, P 2 = x + 1, 

X + 1 
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P 3 = x 2 + llx + 48, P 4 = x 3 + 21x 2 + 218x + 1248, 



(52) 



4.4 Limit (d): The Modular form 

Finally there seems to be one more non-trivial (and, perhaps, attractive) 
possibility. It comes to mind when one looks at the formula (|20|) for Z^,. 
Clearly there should exist a large N limit, when 

z °-ifa (53) 

Naively s = N 2 /2 + 1/24, but one can think about reinterpretation of the 
[/(iV)-dimension N 2 as, say, 

^ +2 ^.£ 1+ /£„^_i +2 (M (54) 

Z n=l n=l Z ^ ' 

i.e. as being regularized with the (^-function technique, thus giving s in 
( |53|) some finite value (s = —1/8). More interesting, the elliptic function 
77(g) = g 1//24 / n^Li(l — Q n ) m the denominator in ( |53| ) has an essential 
singularity as q — > 1: 

I'l-qdexV—!—. (55) 



v(ql) v 1 - qL 

This gives rise to expressions for the correlators like 

k 



which again obey the factorization property in the leading approximation. 
The role of N of the previous cases is now played by N e ff ~ Y^ql = 7E> 
which is large near the critical point; in this sense this limit can be just some 
particular case of (c). 

Note also, that 77(g) in the denominator of (|53| ) is equal to the square root 
of Laplace operator on two-dimensional surface (with two real dimensions, 
q plays the role of the modular parameter, describing the shape (complex 
structure, to be exact) of the torus). 
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4.5 Alternative Large— N Limits Beyond d = 1 



We saw in the previous Sections that the critical behaviour in a Gaussian 
model arose from the presence of zero modes of the quadratic form in the 
action. We also saw that in the Gaussian d=l model exactly the N (= 
rank U(N)) of the N 2 (= dim U(N)) zero modes of the $-field survive after 
the [/-integration. Here, we shall argue that the case of d > 1 is much more 
subtle. 

Let us first consider an analogue of the model C from Section 3 in d 
dimensions. The partition function is given by 



The result of the integral is of course the determinant of the lattice Laplace 
operator raised to the power —N 2 /2. Zero modes on a lattice without bound- 
aries (such as the <i-dimensional torus) are constant (x-independent) matrices 
$. There are thus N 2 of them. On a rectangular lattice they arise when the 
mass vanishes i.e. when the parameter m? = m 2 rit = d. 

If we proceed to the KMM ([!]) the theory becomes non-linear and the 
number of zero modes is no longer N 2 . In fact when d>l it is no longer trivial 
to show that there are N zero modes associated with the diagonal matrices 
<3> (i.e. those belonging to the Cartan subalgebra of U(N)). This is already 
clear because diagonal matrices are not a priori distinguished in eq.(^). In 
the one-dimensional model D they were distinguished but only after the U- 
matrices were diagonalized. Then the diagonal components of $ decoupled 
from [/-fields (since \Uu\ 2 = 1 for diagonal Z7y), and the N corresponding zero 
modes survived the integration over U. Unfortunately there is no analogue 
of this procedure in the higher-dimensional case. Gauge transformations (0) 
are not enough to make all the f/'s at all the links diagonal. (The most one 
can do is to eliminate all the U's on any maximal tree of the lattice and 
diagonalize one more U.) It is however possible to diagonalize <3> at all the 
sites. (The difference is that for d>l there are more links than sites). As 
a result the [/-integral depends only on the eigenvalues of $. The resulting 
formula O] is given by 




(57) 



/($,#) = Vat 



de%)e^ 



V N 



Ep(-l) P e^L^w) 
A(0)A(V>) 



(58) 
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where A(0) = Ylf^ifa — <f>j) an d the sum is over all the N\ permutations P 
of the N variables {ipi, • • • , iPn}- 

Despite the above argument we still claim that there are N zero modes 
even when d>l. To understand the idea of the derivation we begin with 
a simple model which imitates some of the features of the rf-dimensional 
lattice. We consider the KMM on a lattice consisting of a single site with 
d links attached to it. This "lattice" is sometimes called a "bouquet" and 
is shown in Fig. 2. After the [/-integrations the partition function for the 
KMM on this lattice is 



Z D [bou q uet] = / £R P e^-> . (59) 



(A(0)) 2 ( rf - 1 ) 



This integrand is non-singular as fa — <f) 3 ; — > since the order of the zeros in 
the numerator is more than enough to cancel the zeros in the denominator. 
Zero modes of the $-field contribute to the integral at large fa. 

To find the dominant singularity let us simplify the problem even more 
by considering first the case N—2. In this case the eigenvalues of $ can be 
written as fap — 773 (0tr ± V 9 )- The variable fa T which is associated with the 
U(l) factor of U(N) decouples (for any N) and gives rise (for any N) to a 
factor 



/ 



dfa T e- (m - d) ^ ~ ; (60) 
^N{m 2 - d) 



in the partition function. (This is why we usually have one trivial zero mode 
if we consider U(N) instead of the SU(N). For iV=2 the contribution of tp 
to the partition function is 



Zd [bouquet] ~ J d(pe 171 ^ 



_ 2 (e^ 2 -e~^ d 



«p2(d-l) 



The first singularity which arises as m decreases from infinity is at m = 
m^ rit = d. Note that, unlike the situation in d — 1 where the partition func- 
tion itself is singular at the critical point, here it is non-analytic so that high 
enough orders of its derivatives are singular. In d=4 the only the third and 
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higher order derivatives by m 2 become singular, the standard behaviour of a 
third order phase transition. This non-analyticity comes from the large ip be- 
haviour of the integrand which contains the factor e^ d ~ m ^ . It is associated 
with the limit where large <fi is damped only by a logarithmic potential in the 
action, rather than a Gaussian. This is a result of the y? 2 ( d_1 ) term in the de- 
nominator which comes from the Van-der-Monde determinants. We clearly 
see how the [/-integration smooths the singularity at m 2 = d: If this were a 
Gaussian Hermitean matrix model the Van-der-Monde determinant would 
appear in the numerator rather than the denominator and the integration 
would really be singular at m 2 = d. 

It is clear that the partition function can be made more singular at m 2 = 
d by inserting a power of ip in the measure for the $ integration. This 
option, which in the general case amounts to inserting a power of det $ in 
the integration measure, is discussed in more detail in Section 5. 

It is straightforward to generalize this to arbitrary N. The main sin- 
gularity comes from the contribution of the identity permutation P = I 
in the sums fl59|) on every link. The corresponding contribution is ~ (m 2 — 
^jv(iv-i)(d-i)/2-]v/2_ ( We have include the factor of U(N)). As a byprod- 
uct this calculation reveals an interesting feature of the model. Despite the 
non-linearity of the Gaussian KMM the critical value of m 2 is just the same 
as in the quasiclassical approximation. This observation is consistent with 
our hopes that this model may be exactly soluble. 

Note that one should not be confused by the fact that the term P = I is 
only one of AM terms in the sum (|59| ) and that there is also a product over 
all of the links of the lattice. In fact the situation is similar to evaluating 
(1 + (AMe^) -1 )* 01 lmks and it may seem that the second term is negligibly 
small. We saw however in the previous section that this does not happen. 
(The factor of N\ is hidden in the U^ii 1 ~ Q n ) ~ ^K 1 - Q) N as q ~ 1). 
Indeed note that N\e N ~ (Ne) N and in the limit (b), when Ne <C 1 this term 
obviously dominates. 

Unfortunately a similar analysis of the most intriguing limit (d) is difficult 
to do, even approximately, in higher dimensions. It almost certainly requires 
an exact solution of the <i-dimensional KMM (which may not be an absolutely 
hopeless problem) and it will be left for future investigations. 
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4.6 The Large Volume Limit 

To end this section we comment on the large L limit and we explain why 
only case (a) of the large N limits considered above has anything to do with 
continuum limit of the KMM. The crucial point is that in this limit one may 
not assume that q L is close to unity (e^ <C 1) whenever q is near 1 (f < 1), 
One can say that the infrared critical point q — 1 is exponentially unstable 
under renormalization group evolution with increasing L. (Recall |J that it 

is qL = M\ ± \J — 1 rather than the mass Mi that is renormalized in 
the simple manner qi = q L .) It thus follows that our considerations in this 
section are not relevant for the large L limit of the theory. This is made even 
clearer by considering the free energy 

N 2 JL 

-\ og Z D = -— log q L + log(l - ql) (62) 

1 n=l 

which in the normal thermodynamical limit should be proportional to the 
size L of the system as L — > oo. This is certainly true for the first term 
on the r.h.s. of the above equation since logg^ = Llogg for any q). But 
log(l — g£) — > as L — ► oo unless q is fine tuned with exponential accuracy so 
that 1 — q'1, which normally ~ (1 — e~ LNe ) can be replaced by nLe <C 1. This 
becomes even more spectacular in the limit (d) where the main contribution 
to the free energy is 

7T 2 /6 7T 2 /6 

log?7(g z ) -> — 

1 - q L Lt 

which is inversely proportional to the size L of the system. Thus for the 
study of the large L limit only the first term in (^) is relevant. This is of 
course the one which is important for the limit (a), which was examined in 
[0] and which leads to the semicircular distribution. 

One could think that this argument is important to explain the irrelevance 
of limits (b)-(d) only in one dimension. Indeed, there is a great difference 
between the KMM for d=l and for d>l. The main feature of the KMM in 
dimensions greater than one is the presence of many closed contours on the 
lattice while there is only one such contour in model D in one dimension. If 
there were no closed contours, all the [/-variables could be gauged away by 
a gauge transformation (^). This very reason makes the large L limit of the 
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KMM in one dimension (where only one matrix U survives) very different. 
For the study of KMM in higher dimensions it is more reasonable to look 
at the 1-dimensional KMM with small values of L ~ 1. However, one can 
easily see that this argument does not help the critical phenomena arising 
due to the zero-modes: all of them become irrelevant in infinite volume in 
any dimensions^. 

5 KMM Beyond One Dimension: From the 
KMM to a Modified KMM 

Another simple example of the KMM which can be analized rather easily 
is the model with N = 2. We already used this example in the previous 
section, as well as in our considerations of the correlators in ||. We shall 
now analyze the mean-field approximation in the <3>-sector of this model and 
see that it is, first, in agreement with the exact solution at d — 1, second, it 
does not change much when d is increased to be greater than one, and, third, 
implies the existence of the nice mean-field in d > 1 KMM. Of course these 
results do not need to mean too much for the large- N limit of the KMM at 
d > 1, and the conclusion of ref.|7] still can be unavoidable in that situation. 
However, we do not see any obvious reason for this in our simplified analysis. 

The basic formula for the study of the <3>-sector of KMM for N = 2 is 
the integral (^1|). While we introduced this formula for a specific bouquet 
lattice, it has a more universal meaning: it describes the effective theory of 
a constant master field on a conventional rectangular ^-dimensional lattice. 
The value of the master field is therefore defined from the equation of motion 
for effective potential 

V e s((f) = m 2 ip 2 + (d — 1) log — (ilogsinhyj 2 , (63) 

7 The only possibility for these types of critical points to survive can be some sophis- 
ticated limits (similar to conventional double scaling limit of the gravity- models), when 
L and TV are adjusted to increase in a correlated fashion, so that Ne~ Le = const, or, 
in other words, when N grows exponentially with the decrease of the lattice spacing a: 
logiV ~ I /a. This is extremely unnatural from the point of view of QCD, where ane can 
rather expect that the correlation (if any) is N ~ 1/loga. 
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which look like 



m 2 H 5— = d coth. (pi (64) 

while 

&V(<Po) . J d d-l\ 

o 2 = 4( f0 • ,2 2 4- • 65 

<V Vsmh ^ <^ / 

If m 2 < m 2 rit = d, this potential has a form of Fig. 3a, it usually has a 
minimum, but needs and needs some higher order terms to be stabilized at 
infinity || (puncture line on the picture); while in the case of m 2 > m 2 rit = d 
it looks like Fig. 3b. and has a clear minimum at ip = ip® and does not require 
anything in order to be stable. 

In the case of d = 1 the logarithmic term in (^4j) disappears and one can 

get: 



1 , m 2 + 1 

and 



Z w = c -v^o) dct -i/2 d2v M 

^-'mean field — ° ucl q , q , 



K - 1) log" 2 g±l /Vl ' S m 2 -l lt>0 

in agreement (up to the logarithmic factor) with the m 2 — > 1 limit of exact 
formula (pUD for the partition function. This is remarkable, since there is no 
a priori reason to rely upon applicability of the mean-field approximation 
for N = 2. This opens the possibility to believe that analogous results can 
have something to do with reality at d > 1 as well. 

The explanation of the picture Fig. 3b is in fact very simple and seems 
essentially independent of N. From formula ( (59]) it is clear that for in- 
finitely large and different values of the eigenvalues the potential is grows 
quadratically as long as m 2 > m 2 rit = d. When two eigenvalues come close 
to each other, the potential grows logarithmically and they repel. This is 
because of the obvious singularity, arising from A(0) in the denominator in 
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( |59|) (and could be interpreted as attraction of eigenvalues for d > 1) is in 
fact overcome by the zero (of the order 2c?) in the sum in brackets in (|59D (or 
the numerator in Q51f). It looks very plausible that this behaviour, which is 
seen clearly in the case N = 2, shown in the Fig. 3b (repulsion at coinciding 
(p^s and quadratic well at infinity) should persist for all N. 

However, the properties of this minimum are drastically different at d — 1 
and at d > 1. Indeed, at d — 1 the value of the master field (p tends to 
infinity as one approaches the critical point m 2 lit = d from above, while 
d 2 V ((fo) / dip 2 tends to zero, thus giving rise to singularities and the critical 
behaviour of the partition function, but this is no longer true if d > 1. The 
logarithmic term +(d — 1) log0 2 in the effective potential ( |63"D forces it to 
grow at large values of <p 2 even when m 2 = d. Thus the position of the 
minimum remains finite, the curvature d 2 V(ipo)/d(f 2 at the minimum does 
not vanish even for infinitely small m 2 — d and no critical behaviour occurs. 
This makes the result of ref . |7] very natural. The earlier claim of ref . [[J that 
higher order terms in the ^-potential can be adjusted to produce critical 
behavior means that the coefficient of the tr$ 4 term in the potential can be 
adjusted to make the curvature small. 

This, however, does not seem to be the simplest possibility to overcome 
the problem, found in [^]. Our interpretation of the loss of criticality at 
d > 1 attributes it to a logarithmic term in the effective potential (which 
arises from the Van-der- Monde determinants A(0). This effect can easily be 
compensated by adding a logarithmic term to the bare potential. This is 
equivalent to introducing (det$) a in the measure of integration of original 
KMM (|): 

Z D = J d$(x)[dU(x,y)}(det$(x)) a exp ^-^trK($(x)) + 

+ tr$(x)U(x,y)$(y)rf(x,y)). (68) 

<x,y> J 

Clearly a should be non-negative, a > 0, in order to leave the integral over 
Hermitean matrices $ well defined. This is exactly what necessary in order 
to restore the critical behaviour (for m 2 — > d + and d > 1): we just need a 
to be positive (and big enough). In fact for iV = 2 the shape of the effective 
potential in Fig. 3b remains the same for a ^ 0, only the increase at the 
origin is a bit faster: — \og<p 2 — > —(a + 1) log<ft 2 , while the asymptotics at 
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infinity (d— 1) log0 2 + (m 2 — d)4> 2 — > (d — 1 — a) log0 2 + (m 2 — <f)0 2 , is clearly 
dominated by the 2 term if a > d — 1. For N > 2 the picture becomes 
more complicated, since det$ and A(<3>) are very different functions of the 
eigenvalues of $. (Also note, that the trace part of $ does not decouple as 
a ^ 0.) However, for big enough a there usually will be critical behaviour 
at some minima with large and non-coinciding tp^s. 

Indeed, one can think about the entire set of eigenvalues as of a system 
of N point particles on a line (a variant of the Dyson gas). These particles 
interact through a logarithmic repulsion at short distances and a logarithmic 
attraction at long distances (d — 1) log(0.; — 4>j) 2 : a rather peculiar picture of 
a Coulomb gas with a logarithmic "hard core" which is also Coulomb-like. 
The particles also interact with the origin through a combination of a cen- 



tral harmonic potential and a logarithmic repulsion, (m 2 — d)(j) 2 — |-log0 2 . 
Furthermore, since the matrices are traceless, the center of mass of the parti- 
cles is constrained to be at the origin. (Actually, we could consider matrices 
with a non-zero trace but, as we have discussed previously, the trace is an 
irrelevant degree of freedom.) 

As long as m 2 — d > the longest range part of the interaction is har- 
monic and the particles are confined in the vicinity of the origin with size of 
the order 1/ \/m 2 — d. When m 2 — d = 0, the particles still have a Coulom- 
bic long-ranged attraction and if a is not too large, they still form a stable 
cloud. When a is large enough, the repulsive central potential destabilizes 
the cloud and repels the particles to infinity. This is the appropriate critical 
behavior. This system is a rather simple object both for gedanken experi- 
ments and computer simulations; moreover, as we tried to demonstrate, the 
most interesting qualitative results can be extracted from the study of the 
simplest cases, including N = 2. 

The simple analogue of the model 



J d$(det$) Q e- trV ^ +trA *, (69) 

with a single $-matrix and a slightly different interaction with the matrix- 
valued background field A (analogue of U) , which is known as the Generalized 



Kontsevich Model ||, have been analyzed in ||10|| . At least for this model the 
effect of non-vanishing a is known to create no problems for the solvability 
of the model, and is, moreover, natural from the view of its integrability 
structure. Introduction of the parameter a even proved useful for description 
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of continuum limits of the model ||10|| . All this can easily appear true for the 
modified KMM (|68p, where the dependence on the potential V(Q) should 
be studied at least with the purely theoretical motivations (to reveal the 
intrinsic integrable structure, for example). This, however, is the deal for the 
future. 



6 Conclusion 

This paper was devoted to analysis of the <3>-sector of the "Gaussian" KMM. 
The really interesting properties of the KMM (see 0) are essentially indepen- 
dent of the structure of the $-sector. It is only necessary that it guarantees 
the existence of some master-field <fi and that the magnitude of the master 
field is large. We saw that these requirements were fulfilled in all of the 
situations we discussed above, even if the critical point has nothing to do 
with the continuum limit. 

It is of course important to understand the variety of possible universality 
classes associated with the critical behavior. This is why it is also interesting 
to investigate the $-sector of the model. A priori in the KMM it is not 
forbidden to look for a non — trivial universality class represented by a rather 
simple "Gaussian" model. 

We showed that this model indeed exhibits a rich pattern of critical be- 
haviour even at d = 1, but, as we demonstrated, only the most naive large- N 
limit has anything to do with the continuum limit: other types of the critical 
points (defined as the places in parameter space where the partition function 
is non-analytic) do not survive as the volume of the system become infinite 
(they are infrared unstable critical points). Thus, unconventional large N 
limits cannot save the d > 1 Gaussian KMM from the conclusion of 0. 

Therefore, it is necessary to adopt another line of reasoning. In Section 5 
we proposed a simple explanation of the result of |7[|: in the Gaussian model 
and when d > 1 the effective potential for the eigenvalues of <3> remains 
convex for all N. This was due to a logarithmic attraction of the eigenvalues 
at large distances. We can estimate the magnitude of the eigenvalues at their 
equilibrium positions and it is not so large, probably, of order 1-10, probably 
not large enough for the mean field approximation to be effective for the 
£7-sector of the theory. 

This reasoning also suggests a way to overcome this difficulty We intro- 
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duce a repulsive logarithmic central potential which increases the magnitudes 
of the eigenvalues at their equilibrium positions and, if it is strong enough 
and when the curvature of the Gaussian vanishes, it gives a critical behavior 
where some or all of the eigenvalues take on infinite (or at least arbitrarily 
large) magnitudes. This logarithmic potential can be regarded as a modifi- 
cation of the integration measure in the partition function of the Gaussian 
KMM. This is a rather mild modification of the model and may not affect 
its chances to be exactly solvable. 
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Figure Captions: 

Fig. la. Pictorial representation of the trace of a product of $- 
matrices. Dots correspond to the matrices and labels a±, . . . corre- 
spond to pairs of indices in the matrix product. 

Fig. lb. Using Wick's theorem to compute the correlation func- 
tions we must consider contractions of the matrices, the first kind 
of which we denote by a solid line and is depicted here. This con- 
traction of two matrices identifies the the indices as shown. 

Fig. lc. The other contraction, denoted by a wavy line, identifies 
the indices as depicted here. 

Fig. 2 The boquet. 

Fig. 3a The effective potential for SU(2) when m 2 < m^ rit when 

d > 1. 

Fig. 3b The effective potential for SU(2) when m 2 > m 2 CT - lt when 

d>\. 
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